Reversible and irreversible dynamics of a qubit interacting with a small environment 
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We analyze the dynamics of a system qubit interacting by means a sequence of pairwise collisions 
with an environment consisting of just two qubits. We show that the density operator of the qubits 
approaches a common time averaged equilibrium state, characterized by large fluctuations, only for 
a random sequence of collisions. For a regular sequence of collisions the qubit states of the system 
and of the reservoir undergo instantaneous periodic oscillations and do not relax to a common state. 
Furthermore we show that pure bipartite entanglement is developed only when at least two qubits 
are initially in the same pure state while otherwise also genuine multipartite entanglement builds 
up. 

PACS numbers: 03.65.Yz, 03.67.Mn, 03.67.-a 
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I. INTRODUCTION 

Quantum information is providing new perspectives 
and is opening new issues in the analysis of the dynamics 
of open quantum systems. For instance in [l[ the loss 
of coherence of a system of qubits has been shown to 
be linked to the irreversible flow of information between 
the system and the bath due to a build up of entangle- 
ment between the two. Such model has been extended 
to the case in which couplings between the bath degrees 
of freedom are present 0, U El Q ■ In particular it was 
shown that, due to the monogamy of entanglement Q, 
the presence of intrabath entanglement can inhibit deco- 
herence processes. The irreversible dissipation of a qubit 
has been analyzed in 0, H, 0] in terms of an exchange 
of information by means of repeated collisions between a 
single qubit and a reservoir of an arbitrarily large num- 
ber TV of identical independent qubits all prepared in the 
same state. In most of the existing literature an irre- 
versible dynamics is obtained under the assumption that 
the reservoir has a large number of degrees of freedom. 
For instance, in the derivation of master equations the 
Born-Markov approximation is justified on the basis of a 
weak coupling with a basically unaffected large reservoir. 
On the other hand the dynamics of a system interact- 
ing with a small environment can have rich new features 
ranging from memory effects, Poincare recurrences etc. 
. Our work is somehow in the same spirit of the works 
by Mahler and coworkers [Io| , in which the interaction of 
various quantum systems with finite environments, like 
finite quantum networks or quantum Turing machines, 
has been analyzed with particular focus on the onset of 
quantum chaos. In the present paper we will address the 
problem whether, and under which conditions, an irre- 
versible dynamics can appear when a system interacts 
with very small reservoirs, a situation often encountered 
in the physics of mesoscopic systems [ll| and of quantum 
information processing [12l . 1 1 31 ] and of growing interest. 



To this goal we will use a repeated collision model to an- 
alyze the time evolution of a qubit interacting with the 
smallest non trivial reservoir consisting of just two qubits. 
In spite of the apparent simplicity of the model, the sys- 
tem dynamics shows interesting new features. Due to the 
small size of the environment, the reduced dynamics of 
the system is characterized by large fluctuations. How- 
ever an irreversible dynamics is retrieved when the sys- 
tem dynamics is averaged over a sufficiently large number 
of collisions (as we will discuss later this is different from 
the usual time coarse grained introduced in the deriva- 
tion of a master equation). We will show that a - time 
averaged - equilibrium state is reached only for a random 
sequence order of collisions. This is not a-priori an ob- 
vious result. It is important to stress that the reduced 
dynamics of the system qubit cannot, in this case, be de- 
scribed in terms of a Markovian master equation, as we 
will show below. 



II. THE MODEL 

To set the scenario and to illustrate the approach of our 
work let us review the repeated collision model. Consider 
a set of 7V+1 qubits, the first of which is the system qubit 
and the remaining N are the reservoir. The interaction 
between system and environment is due to pairwise col- 
lisions between the system and a singe reservoir qubit. 
Each collision is described in terms of a unitary operator 
Ui. After t collisions the overall state of the system plus 
reservoir is 
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where qse is the total density operator and the sequence 
of labels i\ ■ ■ ■ it specify the order with which the environ- 
ment qubits collide with the system one. In general each 
collision modifies the entanglement between the system 
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and the environment qubits, in particular initially sep- 
arable qubits will become entangled. Such model was 
considered to analyze the processes of thermalization [3] 
and of homogenization 0, Q. These were shown to be 
achieved when the collisions are described by the partial 
swap operator 

Ui = cos??! + ismrjSi, (2) 

where Si is the swap operator between the system qubit 
and the ith environment qubit defined by the relation 
Si\tfj) ® \cf)) — \4>) ® \ip) (here \tp) £g> \(j>) is an arbitrary 
product of single-qubit pure states). The environment 
was assumed to consist of a set of N qubits all initially 
prepared in the (generally mixed) state £ . For instance in 
[3] £ is a thermal state of a single spin. The initial density 
operator was therefore g^ E = g^ The number 

N of environment qubits was assumed to be infinite in 
order to neglect repeated collisions between the system 
qubit and the same environment qubit. Under these as- 
sumption the system qubit was shown to relax to ther- 
mal equilibrium. More in general such model describes 
a homogenization process in which the qubit system ap- 
proaches state £ [8[. Such relaxation mechanism can be 
understood easily if one notes that this models contains 
all the ingredients used in standard derivations of a mas- 
ter equation: the system interacts with an environment 
which rapidly " forgets" the effects of the coupling with 
the system and returns to its stationary state. 

Here we will look at the model from an entirely dif- 
ferent viewpoint analyzing the system dynamics in the 
opposite limit of the smallest non trivial environment, 
namely the limit N = 2. The reason of this choice is 
twofold. On the one hand, although thermalization or 
homogenization are clearly not possible for such small 
environment we would like to know if, and under which 
conditions, the system and the environment reach some 
form of equilibrium. Furthermore, since the entangle- 
ment dynamics among the three qubits can be followed 
and characterized in detail, it is possible to analyze un- 
der which conditions bipartite entanglement and genuine 
tripartite entanglement builds up. 

III. SMALL ENVIRONMENT 

The major new feature of the small reservoir limit is 
clearly the fact that the system qubit collides repeatedly 
with the same environment qubits. In our case Eq. ([1]) 
reduces to 

oil = u it --- u l2 u n e ^ u}jJl---ul, (3) 

where the label refers to the system qubit while 1, 2 
refer to the environment qubits. We will not restrict our- 
selves to the case in which g^ 2 = g$® gi® g2 but we will 
consider also the case in which some entanglement is ini- 
tially present between the environment qubits. The order 
with which the t collisions take place is specified by the 
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FIG. 1: Time evolution of the z-component of the system 
qubit Bloch vector i>o, f° r a random (left) or an ordered 
(right) sequence of collisions. The initial condition for the 
environment is the separable state 1 00} (above) or the Bell 
state -^(|00) + |11)) (below). Parameter values: r\ = r-jTT, 

0=§7T. 

string of indices i± ■ ■ ■ i t . We have considered two limit- 
ing cases: a completely random sequence, corresponding 
to the situation in which the system qubit collides with 
equal probability with each of the environment qubits, 
and the regular periodic sequence 1212 • • • 12 • • • , corre- 
sponding to the situation in which the system qubits col- 
lides alternatively with the two environment qubits. 

A. Approach to equilibrium 

In order to characterize the approach to equilibrium 
we will make use of an important feature of our sys- 
tem: the total Bloch vector, defined as B — 5Zi=o^' 
where hi is ith qubit Bloch vector, is a constant of mo- 
tion. If the state of both system and environment qubits 
approach the same equilibrium reduced density opera- 
tor, after a sufficiently long number of collisions we must 
have hi — ► £>/3. In Fig. [1] we show the time evolution 
of the z-component of the system qubit Bloch vector 
for a random and an ordered sequence of collisions. A 
similar behavior is shown by the x and y components 
(data not shown) . The system qubit is assumed to be in 
the state cos (f ) |0) + sin (|) |1) while the environment 
is assumed to be initially in the state 1 00) or in the Bell 
state ^(|00) + 1 11)). Such choice has been made to illus- 
trate the new features of the dynamics for two different 
instances of environment in a pure state with an equal 
reduced density operator for the two environment qubits 
and with different bipartite initial entanglement. 

The figures clearly show that while for the ordered se- 
quence of collisions the z-component of the Bloch vector 
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exhibits periodic oscillations, strong and rapid fluctua- 
tions appear for a random sequence of collisions, due to 
the extremely small size of our reservoir. This is very 
different from what happens in the large reservoir limit, 
where, for sufficiently weak collisions, i.e. for rj small 
enough, the environment qubits are hardly modified by 
the interaction with the system, while, after a sufficiently 
large number of collisions the density operator of the sys- 
tem qubit becomes monotonically arbitrarily close to £ 

0, |. 

In our system a steady state is however reached if one 
looks at a time averaged dynamics. In order to explore 
the emergence of a time-averaged equilibrium state, we 
have evaluated the components of the vector 

<A&i)(i) ~£ W) -fj= &>(t) - |, (4) 

where (A) (t) denotes the time average of the quantity A 
from time to time t and i labels the qubits. In Fig.[2]the 
length |(A6j Z )| of the z-component of the vector (A6j), 
for the system and the environment qubits is plotted, for 
the same initial conditions and parameter values as in 
Fig. [TJ It is important to note that in the time-averaged 
dynamics one looses knowledge of the exact sequence of 
collisions. If the same time average were done in the 
homogenization process described in all sequence of 
collisions would give origin to the same irreversible dy- 
namics. This is not at all the case in the small reser- 
voir limit. Our results suggest that the approach to a 
time average equilibrium state depends on how random - 

1. e. how compressible in the Kolmogorov sense [3] - the 
string identifying the sequence of collisions is. When a 
regular pattern exists in the collision sequence then the 
time-averaged dynamics of the system and the reservoir 
qubits settle on different values. Note incidentally that 
while in the large N limit, thanks to the fact that the 
system qubit collides with "fresh" reservoir qubits, the 
system dynamics can be analyzed in terms of completely 
positive (CP) maps forming a semigroup leading to an ef- 
fective Lindblad time evolution Q , this is not possible in 
our case, since, due to the repeated collisions among the 
same qubits, the state of reservoir qubits involved in the 
collision changes each time and an irreversible dynamics 
is obtained only after time average. We emphasize that 
our time average is different from the typical coarse grain- 
ing introduced in the derivation of a Markovian master 
equation. Such coarse graining amounts at looking at the 
system dynamics over a timescale longer than the time 
needed to the reservoir to reset to its steady state. In the 
repeated collision model with an infinite reservoir model 
the environment sets to its equilibrium value after each 
individual collision as the system collides with different 
fresh reservoir qubits. 

The straight lines in the left plots of Fig. [5] show a 
decay Abi z oc t~ x l 2 for a random sequence of collisions. 
This implies that the cumulative sum obtained by adding 
the vectors Abi(t') = i)j(t') — ^ from time t' = to time 
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FIG. 2: Length |(Abi z )| of the ^-component of the vector 
(A6i), for the system qubit (i — 0, solid curves) and the 
environment qubits (i = 1, 2, dashed and dot-dashed curves, 
respectively). The straight lines show a decay oc t -1 ^ 2 . Same 
initial conditions and parameter values as in Fig. [1] 

t' = t grows oc i 1 / 2 . Such diffusive growth is the same 
as for Brownian motion or for coin-tossing sequences and 
suggests that, given t\ and t 2 , the vectors Abi(ti) and 
Abifo) are in practice uncorrclatcd, provided \t± — t 2 \ is 
sufficiently large. 

Further insight into the approach to equilibrium is 
gained by evaluating the self correlation defined as 

1 T 

C«(t) = lim — — V Ab ik {t')Ab ik (t' + t), (5) 

T— >oo 1 -j- 1 ^ — ' 

t'=0 

where k — x,y, z. In Fig. [3J we show the decay of the 
correlation function Ci Z (left plots, corresponding to ran- 
dom sequence of collisions) versus the oscillations of the 
same correlation function in the case of a regular se- 
quence of collisions (right plots). The insets show the 
exponential decay of \Ci Z | when the sequence of collisions 
is random. Note that \Ci Z \ eventually oscillates around 
a value oc T" 1 ' 2 due to the finite sequence of data con- 
sidered in computing the correlation function (in these 
plots, T — 10 6 ). This is again compatible with normal 
diffusive dynamics. 

B. Entanglement dynamics 

Interesting features of our model appear when the time 
evolution of bipartite and multipartite entanglement is 
analyzed. To this end let us first introduce the tangle 
as a measure of bipartite entanglement. Given the den- 
sity operator p a b of a bipartite system of two qubits, the 
tangle r a | fc is defined as 

T a \b(p) = [max{0, ai - a 2 - a 3 - a 4 }] 2 , (6) 
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FIG. 3: Decay of the correlation function d z (left plots, corre- 
sponding to random sequence of collisions) versus oscillations 
of the same correlation function in the case of a regular se- 
quence of collisions (right plots). Same meaning of the curves 
as in Fig. [2] initial conditions and parameter values as in 
Fig. [U Insets: |Cj z | vs. time (semilogarithmic plot). 



where {cti} (i = 1, ..,4) are the square roots of the eigen- 
values (in non-increasing order) of the non-Hermitian op- 
erator p = p(o y ® &y)p*(o'y ® <J y ), <7 y is the y-Pauli oper- 
ator and p* is the complex conjugate of p, in the eigen- 
basis of cr z ® a z operator. In our model the tangle Tj\k 
can be used to quantify the entanglement between the 
pair of qubits i,j for an arbitrary reduced density oper- 
ator pij. Furthermore, when the overall state of the sys- 
tem is pure, the amount of entanglement between qubit 
j and all the remaining can be quantified by the tangle 
T i|rest = 4detpj. We have numerically computed the 
tangles r |i, t |2, and T\\2 of the two-qubit reduced den- 
sity matrices and the three-tangle Tjuifc = Tju/.— Tuj— Tjifc, 
where i, j, k can take values 0, 1, 2 and where the tangle 
r^jk measures the entanglement between the zth qubit 
and the rest of the system i.e. qubits j, k. The three- 
tangle t |i|2 is a measure of the purely tripartite entan- 
glement and is invariant under permutations of the three 
qubits @. Our numerical results of Fig. [4] show that, for 
a random sequence of collisions, the time-averaged tan- 
gles (t |i), (tq|2), and (rm) saturate to the same limiting 
value. This confirms once more the approach to statis- 
tical equilibrium for our system. Moreover, we can see 
from Fig. [4] that genuine multipartite entanglement, i.e. 
( T o|i|2) 7^ 0, is created when the two environment qubits 
are initially prepared in a Bell state, but not if they are 
described by the same separable initial state. More gen- 
erally, purely bipartite entanglement builds up among 
the qubits when at lest two qubits (cither both environ- 
ment qubits or system and one environment qubit) are 
in the same state and the initial state is a pure separa- 
ble state. Indeed, if the initial state I^se) is \ip(jxj)), or 
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FIG. 4: Time-averaged tangles (t m) (solid curves), (r | 2 ) 
(dashed curves), (t^) (dot-dashed curves), and (roim) (gray 
curves). Same initial conditions and parameter values as in 

Fig.m 

l^'i/'^i), or \<fxf)il)), then after f collisions we obtain \iPse) = 
a(t)\ip(j)(j)} + (3{t)\4np(j)) +j(t)\(j)(j)ip}, where the coefficients 
a, /?, 7 are determined by the collision sequence. Using lo- 
cal unitary transformations only i.e. with a tensor prod- 
uct of single qubit unitary transformations, we can map 
this state into a(f) 000} + 6(f) 1 100} + c(f ) |010) + d(t) |001), 
where a,b,c,d> and a 2 + b 2 + c 2 + d 2 = 1. This latter is 
the standard form for the states of the W class [l5[ , which 
are characterized by purely bipartite entanglement. How- 
ever, for general initial conditions, genuine multipartite 
entanglement builds up, as is the case when the two en- 
vironment qubits are initially in a Bell state. We point 
out however that the initial presence of bipartite entan- 
glement in the reservoir is not a necessary condition for 
the appearance of multipartite entanglement. Indeed, as 
should be clear from the symmetry argument outlined 
above, a non zero time averaged three tangle (Toim) 
builds up as long as the three qubits are initially in dif- 
ferent pure states. On the other hand if the initial state 
of the three qubits is of the GHZ family the multipar- 
tite entanglement would remain purely tripartite. This 
of course is due to the fact that the partial swap operator 
does not change the overall number of |0) and |1) states. 

IV. CONCLUSIONS 

In summary, we have shown that relaxation (in time 
average) to statistical equilibrium is possible for a system 
of just three qubits undergoing purely unitary evolution, 
provided that randomness is present in the sequence of 
pairwise collisions. We point out that the results of this 
paper can be easily extended to the case in which the en- 
vironment consists of a finite number N of qubits. In par- 
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ticular, a time-averaged equilibrium state characterized 
by purely bipartite entanglement is approached if the ini- 
tial state is separable and N qubits are in the same state 
(either all N environment qubits or the system and N — 1 
environment qubits). In this case, the N + 1-qubit state 
evolves in a subspace of states of the W class (spanned 
by the "ground state" 1 0) and the "single-qubit ex- 
citations" |10...0), 010...0), |0...01)), whose dimension 
is N + 1, much smaller than the overall dimension 2 N+1 
of the Hilbert space. Therefore, the equilibrium state 
can be attained even though genuine multipartite entan- 
glement is not developed. This implies that in a small, 
unitarily evolving quantum system statistical relaxation 



is possible even though the dynamics is restricted to a 
subspace whose states are characterized by purely bipar- 
tite entanglement. Since randomly chosen pure states in 
a many-qubit Hilbert space typically exhibit large mul- 
tipartite entanglement 16], it follows that relaxation to 
equilibrium is possible without exploiting the full com- 
plexity of the Hilbert space. 
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